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Abstract. I study the effect of scalar and spin-orbit rescattering terms, in the 
production of a nonzero Sivers-like asymmetry in proton-proton collisions (inclusive 
production and Drell-Yan) at moderate center of mass energies y/s < 15 GeV 
and transverse momentum up to 3 GeV/c. An ultrarelativistic generalization of 
the Glauber formalism is here used to (i) fit the scalar and spin-orbit interaction 
terms on proton-proton elastic scattering data, including analyzing power, (ii) 
transfer such information to inclusive proton-proton scattering. It is shown that the 
phenomenological interactions responsible for the nonzero analyzing power in proton- 
proton elastic scattering produce a relevant nonzero analyzing power in inclusive 
processes associated with proton-proton collisions. This could represent a relevant 
(possibly higher-twist) contribution to the Sivers asymmetry. 
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1. Introduction 

1.1. General background 

The problem of the study and measurement of T-odd distributions in hadron-hadron 
scattering has recently acquired a certain relevance and quite a few related experiments 
have been thought or scheduled for the next ten years [H [2j [3j HJ E] . 

In particular several studies and models have been proposed for the Sivers 
distribution function[6j. Its possible existence as a leading-twist distribution was 
demonstrated jH [9j [10] recently, and related [TTJ to previously studied T-odd 
mechanisms |1 21 [13]. Some phenomenological forms for its dependence on x and kx 
have been extracted pH M, M, EH EE] from available data [HI M, [2D 1221 [23] . 

While studies of general properties^ E5lEElEZlESlESlEQll3Ill32ll33ll31lES!ofT- 
odd functions relate these functions with a wide spectrum of phenomena, quantitative 
models mostly follow the general scheme suggested in [7]. A known quark-diquark 
spectator model[36] is extended by including single boson exchange [3?] [38j EU HQ]. 
In the case of [H] the unperturbed starting model was a Bag model, and for [12] a 
constituent quark model. 

Here, I want to follow a different approach, that may be of interest in the case of 
intermediate energy measurements, and that relates Single Spin Asymmetries (SSA) in 
inclusive processes with SSA in exclusive processes. 

A well measured SSA in an exclusive channel is the normal vector analyzing power 
measured in elastic proton-proton scattering at energy 20-30 GeV and Qt = 1—3 GeV/c 
(see 031 Si Sana 113 S81). 

In the energy range 20-28 GeV, and for Qt— near 3 GeV/c, this observable is 
unexpectedly large and roughly energy-independent (see fig. 8 later in this work). At 
energies over 30 GeV or Qt over 3 GeV/c it is not measured. 

A nonzero analyzing power in pp — > pp requires interference between helicity-flip 
and helicity non-flip amplitudes with 90° phase difference (see e.g. [19] or [13] for reviews 
on this and related points). For this reason, in PQCD helicity-conserving processes this 
analyzing power must be very small, and the origin of the phenomenon has not yet a 
commonly accepted explanation. Attempts to explain it may be traced back to ref. [51] 
in a non-QCD context. For the region of semihard Qt = 1—3 GeV/c many nontrivial 
models have been proposed [52j [531 El 

but a "standard" model, in the sense of a commonly accepted model for this effect, does 
not exist at present. It must be observed that for some of the previous models (e.g. [61]) 
the effect should persist at much larger energies than 30 GeV, while other ones select a 
preferential energy range, over which the effect should be suppressed. 

The idea underlying the present work is that the interactions producing such an 
asymmetry are also active in inclusive processes originated by hadron-hadron collisions, 
at the same energy come and transferred momenta. In such a case, they could contribute 
to a nonzero asymmetry of Sivers-like kind. So a scheme can be imagined, allowing for 
information transfer from elastic proton-proton scattering to proton-proton inclusive 
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processes. This is attempted in the following. 

I do not propose nor adopt any model for the proton-proton elastic scattering 
amplitude. But I postulate a general form for a set of amplitudes contributing to this 
process at quark level. These amplitudes become initial state interactions in inclusive 
processes, where their presence admits for a nonzero Sivers asymmetry. 

1.2. This work 

The class of processes I want to consider here is the one of SSA in inclusive collisions 
between an unpolarized proton and a transversely polarized proton. In particular, 
azimuthal asymmetries in inclusive production of mesons or (virtual/real) photons. For 
fixed target experiments, the kinematics of interest is the one with moderate beam 
energy 10-60 GeV, and transverse momentum in the semihard regime 1-3 GeV/c. 

To find nonzero T-odd quantities like the Sivers function one needs rescattering 
interactions between the active quarks and the surrounding particles. In particular, 
imagining that a hadron-hadron inclusive hard process is triggered by a hard collision 
between a quark in the projectile hadron (active quark) and a parton in the target 
hadron, this work is centered on the rescattering interactions between the active quark 
and the target hadron (initial state interactions, ISI). 

First, I introduce a very simple bound state in light-cone coordinate representation 
for the active quark in the projectile proton. For this state, I speak of "unperturbed" 
bound state, where "unperturbed" means that it is not affected by ISI with the target 
hadron. This state will be a starting point for reproducing both inclusive and exclusive 
processes. It is a two-component state (each component representing positive or negative 
transverse quark spin). In Appendix A it is shown that in ultrarelativistic conditions 
these two components contain all the relevant independent information on the quark 
state. 

Before the hard scattering event, this state is affected by ISI. The precise form and 
the parameter values for ISI are extracted from the phenomenology of elastic proton- 
proton scattering in semihard conditions. I assume that the interactions producing a 
nonzero analyzing power in proton-proton elastic scattering at the required kinematics 
may be rewritten in terms of interactions between a projectile quark and a target hadron, 
where the former is bound to a projectile hadron, and the latter has a continuous 
structure. 

The scheme used for this aim is an ultrarelativistic generalization of the Glauber 
method[69j. Starting from the fact that the fitted data include unpolarized scattering 
and single normal spin analyzing power, consideration of the number of constrained 
amplitudes reduces the considered ISI to a sum of scalar and spin-orbit interaction 
terms, in a 2x2 formalism. 

Fitting parameters on proton-proton data does not allow for a strict flavor 
separation. However, it implies u-quark dominance. So in the following the expression 
"quark interactions" mainly means "u-quark interactions" . 
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In Section II, some general definitions are presented and discussed, in particular the 
definitions of unpolarized quark distribution and Sivers-like analyzing power, in terms 
of empiric variables on the one side and of quark operators on the other side. Also, 
the quark unperturbed bound state in coordinate representation is introduced. Some 
related details are put into Appendix A section. 

In Section III, an operator describing ISI is introduced. Having to modify the above 
two-component state, this operator consists of a 2x2 matrix operator, that is written in 
eikonal form as the exponential of a set of scalar plus spin-orbit matrices. 

I describe and discuss the theoretical foundations and several details of the method 
allowing me to extract the form and parameters of the ISI operator from elastic proton- 
proton scattering and to apply it to the calculation of inclusive quark distributions. 
Some details are put into Appendix B Section. 

In Section IV the parameters of the ISI operator are fitted to reproduce data on 
proton-proton elastic scattering at beam energy 20—50 GeV and transferred momentum 
1—3 GeV/c, and MRST[70] unpolarized u-quark collinear distribution at Q 2 = 16 GeV 2 . 

In Section V the rescattering operator is used to calculate the unpolarized quark 
distribution, including Sivers-like asymmetry, for x in the valence region and transverse 
momenta up to 3 GeV/c. 

It must be remarked that this work is subject to some limits: 

1) Computational limits. I face stability problems in calculating Fourier transforms 
for transverse momenta over 3 GeV/c (elastic scattering) or 2.7 GeV/c (inclusive at x 
= 0.3). These problems get worse at increasing x, so my analysis centers at x = 0.3. 
This value guarantees that we are in the valence region, and that numerical results are 
reliable. 

2) Limitedness of the data set used for fixing the parameter values. Proton-proton 
scattering is not the only exclusive process from which information on rescattering 
in hadron-hadron hard processes may be extracted, although it represents the most 
complete and precise set of available data in this respect. 

3) Twist identification. The formal apparatus introduced here implies leading twist 
results at the condition that the key parameters of the interaction operator become 
energy-independent at large energies. Since there is presently no way to decide how 
these parameters behave asymptotically, it is not possible to establish whether the found 
effect is a leading or a higher twist one. If it is a relevant higher twist, it should become 
negligible (compared to leading twist effects) at energies like 100 GeV. Exploiting the 
fact that the measured analyzing power is energy-independent from 20 to 30 GeV, I 
assume that the considered interactions are relevant in the energy range 10-60 GeV. For 
what happens over this range, I cannot formulate hypotheses. 

Here the terms "Sivers asymmetry" and "Sivers effect" are preferred to "Sivers 
function". The last one is appropriate in the case of a leading twist contribution. To 
conform with common notation, in the result section I will name "Sivers function" a 
quantity that may be extracted from single spin asymmetries in inclusive processes. 
This quantity must be meant as a measured function, whose theoretical interpretation 
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may be and may be not the one of a Sivers function. 

As a last remark, data and distributions shown in figs. 6, 7, and 8, together with the 
original references, have been reconstructed thanks to the Durham HEP database [71]. 



2. The general scheme - no rescattering. 

2.1. Basic variables. 

Apart for some points where it is clearly specified, all variables will refer to the center 
of mass of the colliding hadrons. Let b = (b x ,b y ) be the quark impact parameter and 
hr = (k x , k y ) the transverse momentum conjugated with it. Let P + be the large light- 
cone component of the hadron momentum, so that xP + is the quark (+) momentum 
conjugated with z_. Assuming by default the validity of a standard factorization 
scheme [721 173] . the fourth coordinate z + plays no role and is fixed to zero. 
I substitute z_ with the rescaled coordinate 

£ = P+z, — ► P + J dz-exp(—ixP + Z-) = J d^exp(-ix^) (1) 

not to work with a singularity of the Fourier transform in the infinite momentum limit 
P + -> oo. 

Contrary to the ordinary treatment of the problem, where one works on a two-point 
correlation operator deriving from a set of squared one-point amplitudes, I develop most 
of the work at the level of one-point amplitudes, square them and then sum over the 
relevant states. 

Since the inclusive process is described here in terms of squared amplitudes, and 
these amplitudes are calculated before being squared, £ is not bound to be positive, as 
it happens in the ordinary treatment based on a two-point correlator with intermediate 
real states. In that case £ has the meaning of the difference between the light-cone 
positions of two points. In this work it describes the light-cone position of one of the 
two only. 



2.2. Two-component transverse spin formalism. 

I consider a quark inside a hadron with a given spin projection S y = +1/2. The quark 
is supposed to present a nonzero transverse momentum along the x direction, and is 
described by a two-component quark spinor 

$ = (M-). (2) 

These two components are the components of the full 4-spinor describing a free 
quark on two 4-spinors ^t+, ^t-'- 

$ = 1p + ^ T + + Ip-^T- (3) 

where 

w {$+,$-), *r- ~ (?-,?+) (4) 
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and 4>± are 2-component eigenstates of 

a y (f± = ±<p±. (5) 

where the scattering plane is formed by the directions "z" (longitudinal) and "x" 
(transverse). 

In the Appendix A section it is demonstrated (i) that in the ultrarelativistic limit 
the above transverse spinors give a full description of the state of a free quark, (ii) that 

^ 7+ * « l^+l 2 + |V>-| 2 , (6) 

that becomes an equality in the u.r. limit. The 7 + — projection produces the distribution 
functions associated to an unpolarized quark with positive z— direction in an infinite 
momentum frame. 



2.3. TMD quark distribution and Sivers-like analyzing power 

I consider the transverse momentum dependent (TMD) quark distribution q(x, h?) of 
unpolarized quarks in a hadron with y— polarization +1/2. This may be defined in a 
phenomenological and in a theoretical way. 

The phenomenological definition adopted here is in agreement with the so-called 
"Trento convention" p3J, in which the unpolarized quark transverse momentum 
dependence distribution has the form 

q(x, k T ) = qu(x, k T ) - ^ - Qs{x, k T ) = qu(x, k T ) ± jjQs(x, k T ){7) 

for an unpolarized quark in a hadron with S = ±y, moving along the z direction. 

If the second term is scale- independent, q$ is the Sivers function. In the following 
I will speak of "Sivers asymmetry" meaning qsk x j1Mqu . 

This asymmetry can of course be isolated by calculating the ratio [q + {k x ) — 
q-(k x )]/[q + (k x ) + q-(k x )] corresponding to opposite proton polarizations. Alternatively, 
one may use the difference q + (+k x ) — q + (—k x ) for fixed proton polarization. The result 
is the same. 

So, the center of this work will be the quark distribution asymmetry A(x,k x ), 
defined as 

\q{x,k x ) + q(x,-k x )J ky=0 

In the following, this inclusive analyzing power of the quark distribution will be simply 
indicated as "Sivers asymmetry", or simply "asymmetry". 

2.4- Sivers-like asymmetry in terms of two-component quark states 

On the theoretical side, q(x, kx) may be defined as the 7 + — projection of the two-point 
correlation function C (x, kx) 

q(x } k T ) = Tr[ 1+ C(x,k T )\, (9) 
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C(x,k T ) = J d£d 2 b e -^k T -b < p|tf( )f(0 + ,f,&)|P > (10) 
oc Y,\M^k T )\ 2 . (11) 

n 

In this work, the strategy will not be a direct calculation of C(x, kx) from the 
definition eqJTOl but a calculation of a restricted number of amplitudes ipn{x,kT) in 
egim These are later squared and summed to obtain the correlator. This choice is 
associated with the need of introducing ISI that have continuous and nonperturbative 
character. 

— # 

The functions ip n = < n\q(x, kx)\hadron > must be read as quark distribution 
amplitudes, i.e. amplitudes for removing a quark with quantum numbers x, kr from 
the initial hadron leaving the spectator in the state \n >J| In any single-particle model 
for a set of bound quark wavefunctions, the above i/j n (x,£) coincides with one of the 
single quark bound state wavefunctions in the infinite momentum frame. The joint 
action of the 7+ —projection, and of the light-cone limit condition z + = (equivalent to 
integration of the quark wavefunction over kJ) selects the correct subspace in the u.r. 
limit. 

If the states ip n of eq|TT] are written in the form of eqJHJ exploiting eqE] I get 

q(x,kr) = Tr[ 1+ C(x,k T )] oc (12) 

oc £Tr(«n7+)= £(V + Mt)| 2 + \i;_(x,k T )\ 2 ) . (13) 



n 



where ip±(x,kx) are complex functions (and not 4-spinors anymore) representing the 
projections of the quark wavefunction on transverse spin states. Inserting eq{13] in eqJH] 
the asymmetry may be written as 

E„(Wx,MI 2 - \4>-{xM)\ 2 ) 

A(x,k x ) = ) (14) 

En(\Mx,kT)\ 2 + |^-(^t)| 2 



2.5. The undistorted quark state 

In this work one state n only is considered, so that the spinor ip = ip-) represents the 
splitting of the (polarized) hadron into a quark with spin projection ±y and spectator 
in this unique state |n >. 

To later insert initial state interactions, we need to express the quark state in 
space-time representation: 



${x,Zr) = J d^d 2 be-^ rb ^,b), ^,b) = (^ + {^) (15) 

I The definition of q(x, kx) in terms of the correlation amplitude is present in several works since (72] 
at least. For more details on its translation in terms of J2 n \tp n \ 2 see e.g. [75] or [76] . 
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I assume that the state 4>(£,,b) is defined by a given value of the the quark total 
angular momentum J in the hadron rest frame, and that in this state a nonzero 
correlation < Shadron • J quark > = +1/4 is present. In other words, the quark J coincides 
with the parent hadron spin. 

A nonzero correlation between the hadron spin and the quark total angular 
momentum is necessary, since a spin-related effect is impossible if a quark transports 
no information on the parent hadron spin. In absence of ISI this correlation would 
produce a nonzero transversity, but not a single-spin asymmetry of naive Time-odd 
origin because of global invariance rules. 

In absence of ISI, we may assume that we are able to calculate the Fourier transform 
eq. f|T5|) and write it directly in impact parameter representation as (PW means "plane 
wave", i.e. ISI-undistorted state) 

iJ(£,b) PW = cf>(0<f>'(b)-\J y = +1/2 > (16) 

A state with given J y may be realized via S and P waves, corresponding to the scalar 
and axial vector spectators of [36J. 

In this work I have avoided, as much as possible, the introduction of parameters 
that cannot be constrained by data. Since elastic proton-proton data may be reasonably 
fitted using a pure S-wave state (see section IV), I have limited myself to the state: 

j&b) PW ee mm ( \ (17) 

where 0(£) and 4>'{b) have Gaussian shapes. 

The width A£ of 0(£) has been chosen comparing the kj- = distribution with the 
shape of the collinear function u(x) as given by MRST[70] at the scale Q 2 = 16 GeV 2 . 

After an initial fit performed in absence of ISI, the fit has been re-tuned again after 
ISI had been included. The final fit is shown in fig. 6, and has been performed with 
gaussian width parameter A£ = 4.5, plus all the later discussed parameters. 

Although the full procedure is recursive, and so no parameter is fully independent 
from the other ones, I may say that the fit on the MRST collinear distribution is 
decisive to establish the width of </>(£), within small corrections. As discussed in section 
IV, longitudinal parameters cannot be extracted from elastic data (this is evident from 
eq J2^|) . On the contrary, the width of <p'(b) and all the rescattering parameters are 
constrained by elastic data (see section IV). 

3. Insertion of initial state interactions (ISI) 

3.1. Basic assumptions of this work 
Assumption 1) 

a set of independent quark-quark interactions is responsible for hadron-hadron 
scattering at E = 10-60 GeV, transferred momenta Qt = 1-3 GeV/c. 
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Assumption 2) 

the same interactions produce initial state distortions of the wavefunction of a 
projectile quark passing through or near a target hadron in hadron-hadron inclusive 
processes, at the same energy. 

Here the word "quark" may also mean "antiquark". Individual quark-quark 
interactions are imagined as exchanges of multigluon sets, that is not trivial to reproduce 
via resummed perturbative calculations. Assumption (1) receives support from e.g. 
the data of [77] for small— Qt scattering, and should be reasonable for Qt ~ 2—3 
GeV/c. Assumption (2) becomes reasonable if interpreted in non-exhaustive sense: 
interactions deduced from elastic proton-proton scattering data constitute one of the 
relevant contributions to ISI in inclusive processes. 

3.2. Basic equations: Scattering between a quark and a composite hadron 

To connect exclusive and inclusive processes, here an ultrarelativistic generalization of 
the Glauber formalism is usedjjjl 

In Appendix B below, a comparison of the steps leading to the nonrelativistic 
and ultrarelativistic forms of the Glauber-distorted quark wavefunction is presented, 
together with a discussion of the most relevant high energy corrections [83j. 

h h h h 




|d 2 b exp[i(k-k')b] 



Figure 1. quark-hadron scattering amplitude in Glauber u.r. approximation. 
The gluon-shaped t-channel exchange reproduces all the terms contributing to the 
scattering, including the no-scattering term that is subtracted on the right side. The 
amplitude for this graph in Glauber approximation (scalar form) is proportional to the 
one in the figure (see text and eq[18|). 

The connection between the amplitudes for describing quark-hadron interaction, 
hadron-hadron interactions, and rescattering distortions of the quark distribution 

§ This set of techniques begins with [69j . For a discussion of the technique, and a review of applications 
to both hadronic and nuclear processes, see [78]. For a detailed example of application to the calculation 
of wave distortions in a non-elastic process, i.e. the case that is closest to what is done here, see [79] , 



'J 



exp i Tfeb) - 1 
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Figure 2. proton-hadron scattering amplitude in the same approximation as for the 
previous fig.l. The main difference, evident in the (scalar) expression for the amplitude 
reported in the figure, is the average over the quark bound state in the projectile 
hadron. For details, see text and ealMl 




d 2 be i(k)b 



exp i 



q/p,n 



Figure 3. The ISl-modified quark distribution amplitude. The corresponding (scalar) 
amplitude is reported in the figure. For more precise details, see text and eg J25I Notice 
that in this case the "—1" subtraction is not needed, since the no-rescattering term 
gives a leading contribution to the amplitude. 



amplitude, is illustrated in figs. 1, 2 and 3. The discussion of this section refers to 
these figures. The (slightly) simplified formulas shown the figures refer to the scalar 
interaction case, while a 2x2 matrix interaction operator is considered in the following 
text and equations. 
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target h 



projectile p 




Figure 4. The full diagram for inclusive production. The blob near the right- 
upper corner is the true hard interaction characterizing the process, also in absence 
of rescattering. The gluon-shaped exchange represents all ISI terms, including the 
no-rescattering one. 




Figure 5. The structure of the ISI, attached to the hard scattering blob. T- 
channel double- line exchanges represent non-reducible sets of gluon, multigluon, meson, 
pomeron, etc exchanges. S-channcl double lines represent intermediate states where the 
projectile quark (or the target hadron) is converted into a more complex intermediate 
state, like q+g or q+qq. In absence of these intermediate excited states the rescattering 
operator would be hermitean. The excited states in the last step (i.e. attached to the 
hard scattering blob) cause a "direct" imaginary part of the mean field, the other ones 
cause an "optical" imaginary part (see text, section III. 3 and Appendix B). 

As a consequence of assumption (1), the elastic scattering amplitude between a free 
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quark and a bound quark in a hadron (fig. 1) can be written as 



T, 



el, quark 



(k ,q) = f(k ) I dV b x}\exp 



iP T(6,0 



1 }Xo, (18) 



(19) 



These relations describe an average of the quark-quark scattering amplitude over the 
target structure. In eqs. [13 and [19j 

-) The coefficient f(k ) is a convention-dependent kinematic factor. 

-) k D z =p q is the initial/final momentum of the free quark, q is the transferred 
momentum, £ is a longitudinal or a light-cone coordinate, P indicates path-ordered 
integration along a constant— b light-cone path. 

-) The 2-spinors Xc Xf assign the initial/final spin of the projectile quark. These 
asymptotic spinors are defined not to contain the spacetime dependence of the quark 
wavef unctions. In other words, the initial/final wavefunctions for the (free) projectile 
quark have the form Xije 1 ^"^^'^ . (these will have to be generalized later, when we 
bound the projectile quark to a projectile hadron). 

-) T(b, £) is a 2x2 matrix in the normal spin space. T'(k ) is the 2x2 matrix of 
forward scattering amplitudes between a projectile quark with momentum k Q and a 
target quark at rest. T" depends on k Q only. The composition of the target proton in 
terms of individual quarks is contained in the single-particle density function p(b, £). 

The target density function p(b, £) derives from an averaging procedure|69j over all 
the spectator degrees of freedom of the target: 



exp 



dl f (1) 



J d2d3....|Vw S et(2,3, ...)\ 2 exp 



dl U{1 



2,1 



3,...) 



(20) 



(21) 



Since here the structure of the target only enters through the averaged operator T(l), 
in the following I work directly with this operator. For this reason it is proper to speak 
of a "mean field" treatment. 

Eq{18]is pictorially represented by fig.l, where the "gluon-shaped" boson exchange 
must not be meant as a single gluon exchange: it represents the full set of possible 
interactions between a single quark and a hadron, approximated by the exponential 
factor in eqfPSl 

In particular, the exponential factor (represented by the first diagram in fig.l) also 
includes the no-interaction term, as evident from the fact that it becomes unity for T 
= 0. For this reason, it is necessary to subtract the diagonal contribution "—1". 



3. 3. Basic equations: distortion of the wavefunction of a free projectile quark. 

The key point is the possibility to rewrite the first term of the right-hand side of eqJTBI 
(i.e. the term that does not contain the "—1" subtraction) in the form 
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f(k ) I d£ I d 2 b ■ e- l{koi -i- b) x^exp 



iP T(b,Z)d£ 



exp 



iP f(6,e)de 



(22) 



so that we may define the concept of "distorted quark wavefunction" , with the 
£— dependent distortion factor 



exp 



iP T(6,£)d£ 



(23) 



The connection between eqs{T8] and [23] is a standard step in the Glauber 
approximation. Within a nonrelativistic formalism it was first used in [69], where also 
the spin-orbit case was discussed. 



3-4- Basic equations: hadron-hadron scattering 

Now two relevant generalizations are possible. First, we may substitute the wavefunction 
of a free projectile quark with the wavefunction of a quark that is bound to a projectile 
proton with momentum pz ± q/2 (fig. 2). In this case the scattering amplitude is (apart 
for an overall kinematic factor) 



T EL (p,q)= Jdtjd% 

- l}^(£,£) (24) 

and it describes the scattering between a projectile proton and a target hadron, for 
q « p. In the following, this amplitude is used to fit a reasonable form for the T 
operator starting from data on angular distribution and normal analyzing power in 
elastic proton-proton scattering at beam energies 20-50 GeV. 

The role of the "—1" subtraction must be stressed. It means that the leading 
scattering term is part of the interaction eikonal operator of eqf24l In the limit of no 
interaction the above amplitude is zero. If the "—1" factor is removed, in the limit of 
no interaction the above equation gives the nonperturbative expression of the projectile 
form factor < ip(jp + q/2)\ip(p — q/2) >. In this case the eikonal operator may only 
contain second order corrections to the form factor, and the "hard" features of the 
scattering event, if present, must be inserted into the large-momentum tail of the bound 
state wavefunction. In this work we adopt the "—1" subtraction, suitable to describe 
semihard elastic scattering with q « p. This means to assume that the bound state 
itself does not contain hard momentum tails. 



iP + (£,b)e 



iqb 



exp 



iP T(£,b)d£ 
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3.5. Basic equations: DWBA and the modifications of the quark distribution functions. 

The other useful generalization is the application of the distortion factor to the 
calculation of the quark distribution function, illustrated in fig. 3. It exploits the Glauber 
approximation within the Distorted Wave Born Approximation (DWBA) scheme. In 
DWBA a matrix element of the form < f\(A + B)\i >, with A and B interaction 
operators, is written as < f'\B\i' >, where \i' > and \ f > are solutions of the problem 
where the interaction B is excluded while A is considered. This allows for a compared 
study of different processes whenever one may assume that the underlying distorting 
factors A have the same origin. The procedure is of course justified if B is harder than 

aM 

The ISI-distorted quark distribution amplitude is 

Mx,k T ) = [d£ f d 2 be- ix ^ +i ^exp\iP ^ f(£,b)d£ 

J J J —OO 



^ n ,Pw{Z,b) (25) 



where ip n ,pw(£,, b) is the amplitude for extracting a light-cone quark in position 
(z-,z + ,zt) = (-P+6, 0, b) while leaving the spectator in the state |n >, in absence of 
ISI. In this work ip n ,pw(£,, b) coincides with the quark bound state given by eq{T71 For a 
generic nonzero T(£,b), eqJ25l defines the quark distribution amplitude that enters eqs. 
TBI and [UJ in presence of ISI. 

It must be noted that the "—1" factor is missing in eqf25l since the undistorted term 
is a leading contribution in this case. In other words, for T(£, b) = the right hand side 
of eqj25] is anyway nonzero and coincides with the ISI-not-affected quark distribution 
amplitude ip n ,Dw(€, b). 

As in figures 1 and 2, the gluon-shaped factor of fig. 3 represents the eikonalized 
full set of rescattering interactions. Since we are speaking of "rescattering" , these do 
not include the true hard scattering event defining the process class (Drell-Yan, meson 
production, etc). To stress this point, in fig.4 I show the amplitude for the full process 
interesting here. The "strictly hard" event is contained in the upper-right blob from 
which a hard line representing a jet, a meson or a massive photon emerges. The gluon- 
shaped ISI include everything softer than the final hard event. 



3. 6. Complex mean field 

With the aim of arriving to a nonzero Sivers asymmetry, it is important that the 
rescattering operator contains different terms able to introduce phase shifts between 
competing amplitudes for the same process: T = J2T{. A special role in this respect 
is played by anti-hermitean interaction terms. I will simply speak of "real" and 
"imaginary" terms to mean hermitean and anti-hermitean 2x2 interaction matrices Tj. 

In fig. 5 the gluon-shaped ISI of fig.4 are expanded in terms of non-reducible t- 
channel exchanges. Intermediate s-channel states that in fig. 5 are named q', q", etc, are 

|| For a detailed example see [80j . Here, high-energy predictions for electron and proton scattering on a 
nucleus A, i.e. for the processes A(p,p)A, A(e, e'p)A — 1, and A(p, 2p)A, are related within DWBA, and 
the Glauber-style distorted wavefunction is compared with distortions calculated by other methods. 
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in general more complex than quark states (e.g. a quark may split into quark+gluon, 
and the gluon may be absorbed in the next vertex). The presence of imaginary terms 
is related with the additional cuts that may be applied to these intermediate states. 

The imaginary terms may have two origins[82, 69, 83J: direct and optical. 

The direct one is the case where the underlying quark-quark scattering amplitude 
T" is complex, and the process is single-scattering dominated. Since the described 
mean field is an average of T' over the target matter distribution, a complex T' leads 
to a complex mean field. This mechanism is surely present in the problem under 
consideration, and corresponds to the single rescattering graph in fig. 5 when q' ^ q 
or h! 7^ h. 

The optical case requires a two-step transition, and is the case where the initial 
state is regenerated after passing through an intermediate inelastic channel. Then, 
the existence of a cut in the intermediate state produces the imaginary part. This 
corresponds to taking the last contributing diagram in fig. 5 with q" = q, h" = h, but q' 
or h! different from q or h. The effective imaginary part is then present in p(b,£). This 
optical effect may be defined of Gribov's kind[83] (see Appendix B). The well-known 
optical effect of Feshbach's kind [82] implies the formation of bound q + h states but is 
suppressed in the short wavelength regime considered here. 

A complex interaction term implies flux nonconservation. Obviously one expects 
absorption of flux from the initial state. To conserve unitarity, this lost flux will 
be redistributed among all the other accessible states. This implies opposite-phase 
imaginary terms dominating at small and large angles in the elastic channel. At small 
angles, we have some elastic scattering due to diffraction, but the net effect on the 
incoming flux is absorption, meaning a negative imaginary part for the interaction term. 
At large angles, second order transitions "elastic — * inelastic — > elastic" produce elastic 
flux. In a single channel approach, this is reproduced by a complex interaction term 
with positive imaginary part. The signature of this is the minimum in the data of fig. 7 
at Q T = 1.2 GeV/c. 

3. 7. Selected interaction terms 

The determination of T = b)T'(k ) entering eqs. fT8lfT9l is not unique, since different 
T(b,£) may lead to the same T^ib) = f^T(b,£)d£. In addition, this operator derives 
from the averaging procedure eqEUthat is under control for simple few-body systems 
only (see [81] for a borderline example). So, in cases like the one interesting here one is 
obliged to guess and fit directly some model form for T(£, b). 

I assume that the distorting factor F of eq. ff23l) does not present a fast dependence 
on x, so in fourier transforms in may be considered as x— independent. This assumption 
breaks down at small x, but that region is of no interest here since it clearly involves 
a different physics under several points of view. In numerical calculations, the path- 
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ordered exponential operator is approximated by a quasi-continuous product: 

expHP J* f(^',b)dA ~ Y[(l + ifd£) (26) 

where the product starts from a negative and large enough £' value where interactions 
may be neglected, and stops at £. The T matrix is a sum of the kind 

T ^ ] T sca i ar -\- ^ * P spin— orbit- ('^'"0 

where two scalar ant two spin-orbit terms are included, and each term has the form: 



f sca lar,n = $n I J J j P(OPnO)' (28) 



tspin-orbit,n = «« | ^ ^ p(Op'n( b )- ( 29 ) 

All the density functions have been chosen with gaussian form, normalized to p«(0) = 1. 
The longitudinal density is the same for all terms and is equal to the squared quark 
distribution amplitude |0(£)| 2 introduced in eqJTTl In the following 3 transverse density 
gaussian functions p' n (b) will be needed: a soft, a semi-hard and a ultra-hard one. In 
all, I will introduce a soft and a semi-hard scalar term, and a soft and a ultra-hard 
spin-orbit term. The S n coefficients are complex, while the a n have been chosen as pure 
imaginary. See next section for details on the fitting procedure and for their values. 

To see that the 0(a n ) terms are spin-orbit terms, I remark that the above matrices 
act on a basis of eigenstates of S y . This means that S y = <7 3 , S z — <J\, S x — cr 2 (apart 
for a factor 1/2). So we may rewrite: 

Tspin-orut d£ = pp'a(b x S y - b y S x ^P + dz_ (30) 

Since k z oc xP + , we identify the triple product between 6, k, and S in the previous 
equation, and this is equal to L • S. 

4. Fit of scalar and spin-orbit interaction terms. 

For extracting the distortion factor eqj23] from eqfJH applied to proton-proton elastic 
scattering, I have assumed that the projectile hadron participates to both processes 
(SIDIS or elastic scattering) with the same quark "intrinsic" distribution amplitude 
eqUH 

The phenomenology of normal spin observables in proton-proton scattering at 20- 
50 GeV, and Qt = 1—3 GeV/c, is not trivial. Whichever the model, more analogous 
terms must be summed to reproduce it. Indeed, data show interference between at 
least two competing terms in elastic scattering [541 14*3] . and three competing terms in 
single-polarization measurements [441 1451 1461 1471 1 
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Figure 6. Unpolarized quark distribution: MRST[70] collinear distribution xu(x) at 
Q 2 = 16 GeV 2 and model xq(x, kx = 0). The model curve is the one that tends to 
zero at small x. 




Figure 7. Fit of elastic pp scattering data. Fel = dcr/dt. Data come from ref.[84] 



Here I ignore the behavior of the analyzing power in the very soft region (where 
it is nonzero but small), and include two scalar complex terms, and two imaginary 
spin-orbit terms. According to the relations of the previous sections, in all the 
cases any of the interaction 2x2 terms contains an overall scalar factor of the form 
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Figure 8. Fit of the analyzing power in elastic pp data. Ael = [FEh{k x ) — 
F E L(~k x )]/[F EL {k x ) + F EL (-k x )]. Data come from refs. [H EH ES EH EE] ■ 

(Pr + iPi)e%p[— (b / Ab) 2 ]exp[— (£/A£) 2 ], where (/3R + i/3i) is the strength parameter, Ab 
the transverse range and 5£ the longitudinal range. 

Also, the parameters of the quark bound state in the projectile are relevant, since 
this state is convoluted with the target interaction operator in the relevant matrix 
elements. This state too contains a factor of the form exp[— (b / Ab) 2 ]exp[— (£/A£) 2 ], so 
we may speak of transverse radius and longitudinal range for the bound state. 

As evident from eqfJUthe data of figs. 7 and 8 cannot give information on any of the 
longitudinal ranges, of the quark bound state or of the interaction terms. For the quark 
bound state this parameter is constrained by the comparison with collinear distribution 
functions [70] (fig. 6). 

I have taken the longitudinal range of the quark bound state equal to 4.5, and 
all the longitudinal ranges of the interaction terms equal between them and equal to 
4.5/v^. This means that for all densities we have pi(£) = |0(£)| 2 . In other words, the 
same longitudinal range is attributed to all the terms of interest here. 

This longitudinal range has been fixed by the fit in fig. 6 a first time by calculating 
the unpolarized distribution in absence of all rescattering terms. It has been fine-tuned 
a second time after including the leading soft scalar ISI term (see below). It has been 
tuned once again after the remaining ISI terms had been included, but with no effect. 

Data in fig. 7 reproduce pp elastic scattering at beam energy 50 GeV, taken from 
ref.jHl]. The fit curve in this figure shows the left/right-averaged Qt— distribution, 
i.e. for each Qt~ value I report the average of the two values of the distribution 
corresponding to Qt = ±\Qt\- After this average, spin-orbit terms have negligible 
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effect on the fitj^j 

So spin-orbit terms are constrained by the data in fig. 8 only. This figure reports data 
points from several experiments measuring the proton-proton normal spin analyzing 
power at 20-30 GeV beam energy [H HQ M, S3 HE] . 

4-1. The soft scalar term. 

In fig. 7 two regions are evident: a soft and a semihard region, corresponding to Qt 
below and over 1.2 GeV/c. The region below 1 GeV/c is dominated by the soft scalar 
term. 

There is an important constraint on this term. For forward scattering, (i) the ratio 
p of the real to the imaginary part of the amplitude is measured and assumes values 
between —5 and —30 % (decreasing with energy) in the region 20-50 GeV [681 [43j. (ii) 
we are in single scattering regime, and the use of a Born-1 approximation is justified. 
This fixes the ratio of the real to the imaginary part of the related interaction potential. 
Indeed, 

A born , 1 (Q = 0) ex Jv(^b)d£d 2 b (31) 

and 

v(t,b) cx ex P ( / 5 f(e,b)de)-i « / 5 nejw, m 

\ J — oo J J — oo 

So, for the soft scalar term that dominates forward scattering, T(£, b) has the same 
phase as p. 

I have taken Re(T) / Im(T) = —0.2 (averaging a decreasing trend in the range 20-50 
GeV). In practice, this ratio has little influence on the following, I have just assumed 
this value to conform with the known part of the phenomenology. 

Since the central data peak is due to peripheral events, I have assumed that the 
transverse range of the interaction is much larger than the transverse radius of the bound 
quark state. With this assumption, the strength and slope of the central peak fix at 
once the strength of the scalar soft interaction term and the transverse radius of the 
bound state of the projectile quark. The transverse radius of the bound state is 1 fm. 
The strength of the scalar soft term is reported below, together with the strength of all 
terms. The precise value of the transverse range of the soft term has no effect as far as 
it is much larger than the bound state radius. I have taken it equal to 4 fm, but with 3 
or 7 fm things do not change. 

More in general, for any interaction term the (transverse) integral of eqEHis cut off 
by the shorter between the bound state transverse and the interaction transverse range. 
When there is a relevant difference between the magnitude of the two, the precise value 
of the larger is not relevant. It must however be remarked that this property does not 

% They have small but nonzero effect. Indeed, the left-right asymmetry is due to the interference 
between spin-orbit and scalar terms. Squared spin-orbit terms produce Qt— even contributions, that 
in the present case are small for Qt up to 3 GeV/c. 
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transfer automatically to eqf2ol because of the absence of the "—1" subtraction in that 
case. This is clearly a source of ambiguity. In the case of the soft term it is anyway licit 
to guess that it has little effect on the asymmetry of the inclusive distributions, because 
of its reduced range in Qt~ space. 

4-2. The semihard scalar term. 

The parameters of the semihard scalar interaction term are given by fitting unpolarized 
scattering data (fig. 7) in the shoulder region at Qt > 1 GeV/c, and by the need of 
obtaining a nonzero analyzing power at Qt > 1 GeV/c (fig. 8). 

The real and imaginary parts of the semihard scalar term are assumed equal in 
modulus. The sign of the imaginary part corresponds to flux production, and the sign 
of the real part to repulsive interaction. The transverse range of this term is fixed to 
0.45 fm by the shape of the data of fig.7 for Q T > 1 GeV/c. 

A qualitative consideration of this semihard scalar term suggest that it mimics, 
within a single channel formalism, the effect of double scattering terms with intermediate 
formation of large- mass states (see the related discussion in section III. 3 and Appendix 
B). 

For fitting purposes, a nonzero imaginary part is needed, with opposite sign with 
respect to the soft term, to reproduce the interference pattern at 1.2 GeV/c. A nonzero 
positive real part does the same, because of a similar interference with the real part of 
the soft term. Because of the real to imaginary ratio of the soft part, with these signs 
the real part of the semihard potential contributes to 20 percent of the interference effect 
in the dip region. 

In presence of a zero imaginary part, one would require a much larger strength for 
the real part to reproduce the dip. This larger strength would not fit the region Qt > 
1.2 GeV/c. The same would happen for a negative real part: the imaginary part should 
be increased to compensate. 

The major role of the real part is to interfere with the imaginary spin-orbit terms 
producing a nonzero analyzing power. So, the real constraint on the real part arrive 
from the joint fit of data of fig.7 and fig.8. In addition, the strengths of the real part of 
the semihard scalar term and of the hard spin-orbit term are not fully independent. 

4-3. Spin- orbit terms 

Analyzing power data present three characteristic regions, only two of which evident in 
fig.8: (i) a low but nonzero (positive) broad peak at very small Qt (not showed in fig.8, 
but visible e.g. in [44]); (ii) a peak between 1 and 1.8 GeV/c; (iii) a large increase over 
1.8 GeV. 

Here, I neglect the ultra-soft peak at small Qt- I introduce two spin-orbit imaginary 
terms: a soft term, with the same density structure as the scalar soft term, and a hard 
term, with a very small transverse interaction range 0.16 fm. The interference between 
the former and the semihard scalar term (real part) produces the nonzero values in the 
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region 1—1.8 GeV/c, while the interference between the latter and the semihard scalar 
term (real part) produces the large— Qt rise. 

The strength parameter for the soft spin-orbit term has been tuned so to best 
reproduce those points that present the smallest error bars in the confused data set of 
the region 1 — 1.8 GeV/c. 

The hard spin-orbit term is necessary to reproduce the increase of the analyzing 
power in the Qt > 1.8 GeV/c region. In fig. 8, some points at large Qt with large error 
bars have not been reported. If taken at their central value, they would suggest even 
larger analyzing powers than the fitted ones. In addition, my fitting curve decreases for 
Qt > 3 GeV/c. The decrease is driven by the fall of the scalar semihard term interfering 
with the spin-orbit terms. But we do not know what happens to data at Qt > 3 GeV/c. 

The dip at 1.6—1.8 GeV/c is due to the fact that neither of the two introduced 
spin-orbit terms is strong there. I have not introduced a semihard spin-orbit term with 
range 0.5 fm as for the scalar potential, since it would produce large analyzing powers at 
Qt = 1.6—1.8 where we see small or even negative analyzing power. Some of the small- 
error points reported in the figure are born from measurements that were specifically 
dedicated to the region 1.6—1.8 GeV/c. They confirmed with little doubt that in this 
region the analyzing power is close to zero. 

The other dip, at 1 GeV/c, is well measured, and is due to the cancellation between 
the real parts of the two scalar interactions. Then the imaginary spin-orbit soft term 
cannot interfere with anything. On the left of this dip my fitting curve produces a very 
small, negative, rather flat analyzing power that does not correspond to reality (see 
e.g. [15)- Data show an again positive, but smaller than 5 %, analyzing power in the 
most forward region Qt < 0.5 GeV/c. So, an extra term would be needed to explain 
quasi- forward data. Because of the smallness of the effect, I have not cared data at Qt 
< 1 GeV/c. 

4-4- Parameter values 

Bound state longitudinal range: A£ = 4.5. 

Longitudinal range of all the density widths: A£ = 4.5/v / 2- 
The transverse width is peculiar: 
quark bound state: Ab = 1 fm. 

soft scalar and soft spin-orbit interaction: Ab = 4 fm. 

semihard scalar interaction: Ab = 0.45 fm. 

hard spin-orbit interaction: Ab = 0.16 fm. 

The strength parameters are: 

soft scalar interaction: 6 = 0.07 (i-0.2). 

semihard scalar interaction: 5 = 0.01 (1-i). 

soft spin-orbit interaction: a = 0.0004 i. 

hard spin-orbit interaction: a = 0.0001 i. 
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5. Estimate of effects in inclusive processes 

5.1. Distributions and asymmetries at quark level 

Fourier transforms have been calculated assuming kj> coinciding with (k x ,0). 

Fig. 9 shows the quark distribution for positive and negative k x , for x = 0.3. The 
two curves of fig. 9 may also be read as the distributions corresponding, in both cases 
for positive k x , to a parent hadron with spin oriented as ±y. Their non-identity implies 
a Sivers-like asymmetry and a nonzero function qs(x, kr) according with eqE q(x, kr) 
= qu{x, &t) + (k x /M)q s (x, kr), for an unpolarized quark in a hadron with polarization 
along +y. For opposite hadron polarization {k x /M)qs is substituted by —{k x /M)qs (see 
eqI7|). The qs/qu ratio extracted from the difference between the two curves in fig. 9 is 
shown in fig. 10 as a function of kx- 



q(x = 0.3, k x =±k T , k Y =0)/q(0,0,0) 
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Figure 9. k x — dependence for x = 0.3 and k y = of the quark distribution q(x, k x ,0) 
(see eq|7] in the text) for positive and negative k x . These two curves refer to a parent 
proton polarized along +y and to an unpolarized quark. However, they also represent 
the positive k x distribution in the case of parent protons with opposite polarization 
±y. The asymmetry between the shown curves is a Sivers-like asymmetry, from which 
the ratio between the Sivers and the unpolarized distribution may be extracted (next 
figure) . 

Observations: 

1) Figs 9 and 10 are completely representative of what may be seen at other values 
of x in the valence region, e.g. x = 0.2. 

2) Repeated application of the numerical calculation code shows that numerical 
calculations lose progressively reliability for kx > 2.7 GeV/c at x = 0.3. For larger x 
this loss of reliability is more severe, because of the fourier factor exp(ix^). For this 
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Figure 10. Ratio between the Sivers and the unpolarized distribution functions 
qs(x, kx) and qu(x, fcy), as extracted from the asymmetry of the two curves reported 
in the previous figure (see eq[7|in the text). For kx < 1.8 GeV/c the curve does not 
overcome 0.02. 

reason the value x = 0.3 has been chosen, as the largest x value for which a reasonable 
k T — range could be covered. For the been fixed to zero. 

Observing fig. 9, although a flattening of the distribution at hr over 2 GeV/c is to 
be expected, the steep rise of the k T distribution at k T « 3 GeV/c is not reliable. Up 
to 2.7 GeV/c the results of the numerical calculations are stable against changes of the 
integration point numbers and/or integration ranges. In the case of the figures referring 
to elastic scattering, stability arrives to 3 GeV/c. 

3) The peak of the asymmetry evident at 2.4 GeV/c in fig. 10, together with the 
disappearance of the asymmetry over 3 GeV/c, is not an effect of the code unreliability, 
but of the chosen interaction shapes. In other words, although numerical predictions 
are not reliable at 3 GeV/c, the fourier transform of the involved interaction terms must 
decrease in this region. 

4) If one wants to compare the shape of the curves in figs. 9 and 10 with those in 
figs. 7 and 8, one must observe that the transverse momentum scale of any phenomenon 
is slightly reduced in the quark distribution case. This is due to the fact that in eqJ24lwe 
find IVKC; b)\ 2 , while in in eqj25l we have 6). Because of the chosen gaussian shape, 
this means slightly harder effects in the elastic case. 

5) Of the interaction terms inherited from elastic scattering, only the semihard and 
hard ones produces relevant effects on the quark distribution side. Soft scalar and spin- 
orbit terms are necessary to explain elastic data, but their effect on the asymmetry in 
the case of inclusive processes is small. In other words, the peak at 2.4 GeV/c in fig. 10 
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can be reproduced after excluding all interaction terms but the quoted two. For k T < 
1.8 (not shown in fig. 10) the asymmetry is never larger than about 0.2. 

6) The fall of the asymmetry at small kx is an unavoidable consequence of the 
reduced size of the SSA in elastic scattering. For k T > 3 GeV/c it is a consequence of 
the lack of hard interfering contributions, extracted from elastic scattering. It must be 
remarked, however, that in the case of large kx this lack is not constrained by elastic 
scattering data. Elastic scattering data do not show any decrease at large kx- Simply, 
we do not know what happens for kx > 3 GeV/c. But in fitting figs. 7 and 8 I have 
adopted a "minimal" approach, only introducing the strictly necessary interaction terms 
to reproduce the presently available data ranges. Both in the case of elastic and inclusive 
data this approach may lead to underestimation of asymmetries at large kx- 

7) Curves in figs. 8 and 10 present qualitative differences at intermediate kx- This 
is due to two facts: (i) the presence of a fourier transform exp(ix£) in the inclusive 
case, that changes the phase properties of the interaction terms, (ii) eqj25] does not 
contain the "—1" factor of eql24l so in the inclusive case we have interference between 
no-rescattering and rescattering terms. 

8) As specified in the Introduction, it is not possible to know whether the predicted 
effect is a leading twist or a higher twist one. So, if one names "Sivers function" the 
function qs(x, kx) that may be extracted from the asymmetry between the two curves 
of fig.9 using eqd the result shown in fig. 10 may be read as a Sivers function in the 
sense that it mimics a Sivers function at finite energies. 

5.2. Convolution with a partner distribution 

In all the relevant phenomenological cases, the calculated functions must be convoluted 
with a partner function f(x',k' T ), so that the experimentally detected transverse 
momentum is Qx = kx + k' T . 

In a Drell-Yan event / is the momentum distribution of an antiquark or of a gluon. 
In the case of meson production in hadron-hadron scattering, / will be a convolution 
of distribution and fragmentation effects, so that x' and k! T will be product/sum of 
distribution and fragmentation variables. In both cases, we must distinguish the cases 
where k' T derives from Fermi motion [k' T « 0-0.7 GeV/c) or from hard gluon secondary 
radiation (k' T » 0.5 GeV/c). 

The former one will be the case for Qx < 3 GeV/c and for beam energies << 
100 GeV. It is then realistic to imagine a gaussian distribution with pure Fermi motion 
origin like the following one: 

f{x\k' T ) oc e -(*r/*o> a . (33) 

So, fig. 11 in the following refers to this situation. For larger Qx and larger beam energies 
the requirement of a transverse momentum with dominating Fermi motion origin is not 
justified. 

To correctly calculate a convolution for assigned x and x', I need the values of 
q(x, kx) (the function appearing in eqJT] and plotted in fig.9) over a wide 2-dimensional 
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Figure 11. Estimate of a possible phenomenological asymmetry, e.g. in Drell-Yan 
or in meson production. Here the two quark distributions reported in fig. 9 have been 
convoluted with a function f(x', k' T ) = exp[—(k' T /k ) 2 ]f(x') (see text), and Qt = kr + 
k' T . The three curves correspond to k — 0.4, 0.5, 0.6 GeV/c (the asymmetry decreases 
for increasing k a ). The partner function f(x', k' T ) may be the distribution function for 
an antiquark or gluon, or a combination of distribution and fragmentation effects, or 
also include recoil effects by radiation of undetected gluons. As evident, the asymmetry 
is visible if the momentum spread has Fermi motion origin (k a « 1 GcV/c) but not 
in the case where the momentum spread derives from gluon radiation (k D > 1 GeV/c). 




range of kr- As previously observed, the numerical calculation of q(x, kx) is restricted 
by computational problems, that increase at increasing k x and k y . For x = 0.3 and k y 
= 0, I have a reliable set of values of q(x,kr) up to k x = 2.6-2.8 GeV/c. To estimate 
a convolution, I make the simplifying hypothesis that q(x, ±/c T ) depends on k y via a 
factorized term 

q(x,k T ) = q(x,k x ,0)q'(k y ). (34) 

With this hypothesis, the terms that depend on Q x and Q y separate in the 
convolution / dk x dk y q(x, k x , k y )f(x', Q x — k x , Q y — k y ). For Q y = 0, I simply obtain 

/dk x dk y q(x, k x , k y )f(x , Q x — k x , Q y — k y ) — (35) 

= const ■ J dk x q{x, k x , 0)f(x', 0)e~ [iQ *~ k * )/h>]2 (36) 

As for the calculation of fig. 10, the above integral may be calculated for positive 
and negative Q x . The corresponding asymmetry (G+ — G-)/(G+ + G-) is shown in 
fig.ll. 
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Obviously the asymmetry shown in fig. 11 inherits the numerical uncertainties of 
the calculation of q(x, hr), and for this reason it cannot be considered 100 % reliable 
at the largest reported Qt- As for fig. 10, the qualitative origin of the asymmetry peak 
is the interference between semi-hard and hard interaction terms. The semi-hard term 
decreases for kx > 2.5 GeV/c, and so indirectly for Qt > 3 geV/c. So, qualitatively 
fig. 11 shows nothing strange, but it cannot be considered precise for Q T > 2.7 GeV/c. 

Comparing fig. 10 with fig. 11 one notices that in the latter case we have a predictable 
spread of the asymmetry peak towards large Qt, but not towards small Qt- This is 
due to the steep decrease of the quark distribution q(x, k x , 0) (see fig.9) for k x = 2—2.5 
GeV/c. Events with k x = 2.5 GeV/c and k' x = —0.5 GeV/c are much less frequent than 
events with k x = 2 GeV and k' x = 0, so the former kind of events does not influence the 
region Q x rs 2 GeV/c. 

For this reason, if the average fc^— spread of the partner distribution f(x',k T ) is 
increased from 0.4 GeV/c to larger values, the peak asymmetry decreases. So, the 
phenomenological effects are expected to depend in a marked way on the specific 
measurement Q 

Referring to the discussion in point (6) of subsection 5.1, it must be remarked 
that the fact that the found asymmetries decrease for Qt < 2 GeV/c is a necessary 
consequence of the small (average) values of SSA in elastic scattering for Q T < 1.8 
GeV/c (fig. 8). So, in the case of a relevant asymmetry (over 5 %) measured at Qt < 1 
GeV/c, I would exclude that it may have the origin that is described here. 

On the contrary, there are no experimental constraints on the asymmetry values 
from elastic data at Q T > 3 GeV/c, where more interaction terms could be present, 
aside of those considered here. I have adopted a "minimal" approach, in the sense of 
introducing just those interaction terms that are strictly necessary for reproducing the 
visible data. So the prediction reported here could underestimate asymmetries on the 
right side of the peak in figures 10 and 11. 

5.3. Discussion 

As previously observed, I do not imagine the term calculated by me to saturate the 
Sivers function. Rather, it is a contribution to it in a well-defined kinematic region. For 
the following discussion, I will name 'AB-term" this contribution. 

I will discuss the Drell-Yan application, since in this case we have some proved 
factorization statements [721 [73], a simple connection with the SIDIS (leptoproduction) 
case [8], absence of final state effects, a reference hard scale Q (the dilepton mass) for 
Q T , and a connection between the large-Qr and the small-Qx behavior of the Sivers 
asymmetry [llj. 

From figs. 10 and 11 it is evident that the AB-term may be relevant in the region 
A << Q T « Q, for reasonable values of Q. A may be any soft scale parameter. This 

+ In the meson production case things are also complicated by the presence of direct contributions to 
the asymmetry from the fragmentation side [85 j. 
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region has been studied [TT] because here the twist-3 model by refs. [TS1 E3] and the Sivers 
function scheme have overlapping regions of consistency. I name "semihard Qt region" 
the above region. 

The previous calculations refer to proton-proton collisions. As far as it can be 
considered scale-independent, the AB-term should be present in SIDIS too according 
with the change-of-sign rule stated in (8] for the leading-twist Sivers function. 

In Drell-Yan we have two relevant scales: Q and Qt- It is normally admitted [6 J 
that the Sivers asymmetry is power-suppressed in Qt, but the associated Sivers function 
is anyway quoted as "leading twist" if it is not power-suppressed with respect to 1/Q for 
fixed x and A>r- Since Q is related to the squared cm. energy s via the scaling variables 
X\ and x 2 (Q 2 = Xix 2 s)) "leading twist" means that in eqJTJgs depends on s at most 
logarithmically (for fixed x and £>r). 

In this work s is hidden in eqJT9l in the dependence of T' on the momentum k a of the 
projectile quark (in a frame where the target hadron is at rest) This dependence means 
that the strength parameters 5 n and a n in eqs. [28] and [29] are in principle functions of s. 
These parameters are extracted from the data in fig. 7 and 8, so they can be considered 
as s-independent if these data do not depend on the beam energy. The data reported 
in fig. 7 are stable for beam energy 20-50 GeV, and would show logarithmic changes at 
larger energies. Those of fig. 8 are stable in the range 10-30 GeV, but we have no similar 
data at larger energies. So, I cannot presently establish whether the strength of the 
spin-orbit terms is s-independent or is not. If it is, the function shown in fig. 10, with a 
change of sign, is a contribution to the Sivers function in leptoproduction. 

If I decided to calculate directly (with a technique that is reasonably similar to the 
one adopted here) the SIDIS n + production on proton, this would be the crossed process 
of 7r _ -proton Drell-Yan. For the Sivers asymmetry in these two processes the calculation 
performed within my scheme would respect Collins' rule (since the only change between 
the two cases would be in the integration path for the eikonal factor). However, the 
calculation of any of these two processes would require additional assumptions, since 
the elastic data of figs. 7 and 8 only constraint proton-proton ISI. 

An obvious problem with figs. 10, 11 is the non regular rise of the asymmetry 
with respect to Qt- A nonzero AB-term in the semihard region does not exclude the 
simultaneous presence of a Sivers function like those proposed in refs. [371 EHl EHJ H0J 
Bdl W2\ filling the soft Qt region where the AB-term is small. Indeed, we may have 
mechanisms that are able to produce a single spin asymmetry in inclusive scattering but 
become ineffective, or effective but scarcely visible, when applied to elastic scattering. 
These mechanisms would escape the presented analysis. 

On the other side, although nothing in my model forbids a nonzero AB-term at 
small Qt, it is a matter of fact that, whatever mechanism produces a nonzero analyzing 
power in elastic scattering (data in fig.8), this mechanism has its top relevance in the 
region Q T > 2 GeV/c, and rather small effects at Qt < 1 GeV/c. So, to invent a 
model that transforms these small effects, observed in elastic scattering at small Qt, 
into relevant effects in inclusive processes at the same Qt would not be trivial. 
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A over-simplified way to reproduce the physics described in this paper can be 
to imagine that, before the hard inclusive event, the two colliding protons scatter 
elastically and remain almost on shell up to the hard event We know from the data 
of fig. 8 that after the elastic scattering the space distribution of the colliding protons 
depends on the normal spin of the initial state. This is the entrance way to a nonzero 
Sivers asymmetry, since these protons enter the later hard scattering with asymmetric 
momentum distribution. The same data in fig. 8 tell us that this effect is remarkable 
only when the Qt that is exchanged in the elastic scattering is semihard. We note that 
the effect is present also if quarks are completely collinear in the initial state, since what 
matters is the Qt exchanged in ISI. So one has a physical picture where the typical 
event characterized by a nonzero AB-term has small "primordial" Qt, and a semihard 
Qt produced in ISI. 

For the models of the Sivers asymmetry in Drell-Yan that are present in my reference 
list, we may distinguish two classes: 

(A) models where the transverse momentum dependence of the cross section is 
entirely of Fermi motion, "primordial", origin. Here Qt = hri + hri is the sum of 
the transverse momenta intrinsically associated with the distribution functions of the 
colliding quark and antiquark. 

(B) models where the distribution functions are initially collinear, the Qt- 
dependence of the cross section is associated with the hard nucleon-nucleon interaction 
(e.g. gluon radiation accompanying the hard vertex) and some arguments allow one to 
stick (part of) this effect to the individual quark distributions. 

All the models extending the one by [7] belong to group A. The calculation in [TT] 
belongs to group B. The model discussed here is mid- way between the two groups: the 
unperturbed quark distribution is narrow, but not fully collinear (it has a gaussian tail 
and transverse range 1 fm) , and the event numbers at semihard Qt are strongly enhanced 
by ISI. The Glauber-Gribov approximation allows me to stick this enhancement to the 
quark distributions. In absence of ISI the curves in fig. 9 would be lower in magnitude 
for Kt over 1 GeV/c. 

6. Conclusions 

Summarizing, starting from the assumption that the quark total angular momentum is 
dominantly oriented as the parent hadron spin, it is possible to build a nonzero Sivers- 
like asymmetry via mean field initial state interactions of scalar and spin-orbit kind. 

The specific form of the interactions, together with the values of the parameters, 
is here extracted from the phenomenology of proton-proton elastic scattering at 20-50 
GeV. When these interactions are applied to transverse momentum dependent quark 
distributions, they produce an asymmetry that is important at transverse momenta 2- 

* This process is power-suppressed in Qt with respect to the more general process considered in this 
paper, since it requires re-formation of the proton bound state after ISI and before the e.m. hard 
scattering. 
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3 GeV/c, and reduces to very small values for smaller transverse momenta. It is not 
possible to estimate what happens over 3 GeV/c, for the lack of corresponding data in 
elastic proton-proton scattering. 

With the present day knowledge of spin-orbit hadronic interactions, it is not possible 
to know whether their effects persist at very large energies. Consequently it is not 
possible to establish whether a Sivers-like asymmetry generated by them is a leading 
twist one or just an intermediate energy effect. 



7. Appendix A 

Solutions of the Dirac equation have in general four independent components. This is 
related with the existence of two solutions of the scalar equation E 2 = p 2 + m 2 and of 
two independent helicity states for each energy. In the two limiting cases 

2 

E w m + (37) 
2m 

E » Pz + ^, (38) 
2p z 

the structure simplifies. Although the use of a standard 4-component Dirac formalism 
may allow one to exploit a large set of well-known relations for trace calculation etc, 
it hides the fact that two components only are independent in the ultrarelativistic or 
nonrelativistic spinors. 

Writing a 4-spinor in standard representation as 

* = (0,X) (39) 

the two 2-spinors satisfy the coupled equations 

(E - m)4 - p. a X 
[E + m)x — p ■ cr<p 

equivalent to the Dirac equation. 

We have the known nonrelativistic limit \ — > for particles, and the opposite limit 
for antiparticles. This decouples each other spin states, since the second one of eqs 14131 
loses relevance, and in the first one the right-hand side may be neglected. So, only 
interactions may remix spin states. 

As well known, in the opposite limit E » m we have two relevant states </>±x that 
decouple each other. The nature of these states is clear after introducing helicity n ■ a. 
The u.r. limit may be reached for E ~ ±p 2 . For E ~ +p z , <P + X coincides with one of 
the two helicity states, and <p — x with the other one. For E pa —p z , the correspondence 
is the opposite. 
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Instead of using helicity basis, I write the above equations HD] in terms of eigenstates 
of the spin projection along the x— axig^l <j> x ± = 4>t±- 

T _ = (1,1)/ A <h- = (1,-1)/V2 (41) 

These 2-spinors satisfy the overlap relations 

&y4>T+ = ~i<t>T-, CTy4>T- = i<pT+ 

and the same for the corresponding \t± terms. 

In the u.r. limit, from eqH0]one gets the relations 

T+ « XT-i 0t- « Xt+ (43) 
meaning that a spinor polarized along the transverse axis has the approximate form 

= (0 T+ ,XT+), ~ (0T+,0T-), *r- ~ (0T-,0T+) (44) 

More precisely, in presence of a nonzero transverse momentum p y « p z that is also 
orthogonal to the normal spin, eqsfJO] become (exploiting the previous eqsH2]): 

E<p T+ - p z XT- + m<f) T+ + ipyXr- = 0, 

E-XT+ ~ PzXT- + W>T+ + iPyXT- = 0, 

so that 4>t+ and xt- decouple from 4>t- and xt+ (remark: 4>t+ and xt- are pieces of 
different 4-spinors). Writing E = (E—p z ) + p z , the former of the previous two equations 
becomes 

p z (4> T+ - Xt-) + (E-p g )4> T + + ni(j) T+ + ipyXr- = 0. (46) 

At 0(1/ E) this confirms the equality between <p T+ and Xt-- Substituting xt- = 
(pT+, and exploiting E — p z ~ m^/2p — (m + ip y )(m — ip y )/2p one may also quantify 
the difference between the two spinors, that is 0((pT+P y /Pz)- 

For the 7+ projection we have (proportionality factors substitute powers of l/\/2) 

^+ +7 o7+^t+ cc ^+ + ^ J ^ J = (47) 

OC |0 T+ | 2 + <Pt+XT- + XT+0T- + \XT+\ 2 OC |0r+| 2 + |0T-| 2 - (48) 

and the same for ^t-- In addition, 

#£ + 7o7+*t- = #£_7o7+*t+ = 0. (49) 
So, if a generic 4-spinor \1/ has the form 

* = V+^T+ + 1p-^T- (50) 

we have, within corrections 0(1/ E), 

« * + 7o7+^ « IV'+I 2 + l^-l 2 , (51) 
i.e. the result exploited in eq{T3]of the present work. 

(j In this part of the work I project the transverse spin along the a;— axis and the transverse momentum 
along the y— axis. This simplifies intermediate passages if one wants to reproduce oneself the above 
relations. In the rest of the work I project the spin along the y— axis to conform to ordinary treatments. 
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8. Appendix B 

The exploitation of the Glauber technique in DWBA is better suited in nonrelativistic or 
ultrarelativistic form. In both limits it is often possible to approximate the kinetic energy 
in such a way that (i) a large part of the boost energy is removed from the problem, 
(ii) interactions may be inserted as additive contributions, to the zero or to the plus 
component of the energy-momentum vector. When (i) is true, the two treatments are 
equally effective. When it is not so, corrections are needed in both cases[83j. 

In the n.r. case, one normally uses a linearized Schroedinger equation, but 
relativistic kinematics. Since one always works with momenta, and interactions are 
directly fitted within this scheme, the hamiltonian is never involved and it is not 
surprising that this scheme may work in relativistic conditions. A detailed discussion of 
this and related references may be found in [80] . 

Limiting the discussion here to the scalar case for simplicity, in the n.r. case one 
of the possible starting points for the Glauber treatment is a linearized form of the 
Schroedinger equation, via the approximation 

(k 2 -p 2 )^ » 2k(k-p z )?p (52) 

where the kinetic energy is (hk) 2 /2m. Rescaling the potential energy as 2kV = 2mU /% 2 
the linearized Schroedinger equation is 

p z i> = k{l + V)ip, V « 1. (53) 

with stationary solution (for V « 1) 

i) in = exp(ikz + iJ Z V(b,z')dA, (54) 

that is the n.r. version of the wavefunction ip appearing in eqJ22l The corresponding 
solution il) out has Jf^ substituted by . 

For the calculation of the matrix element < ipf\ip > in eq f22l the eikonal direction 
z of the dominating exp(ikz) factor is rotated to kz ± (q/2)x in the initial or final wave, 
while the eikonal distortion factor is still approximately calculated along the direction 
z that is the average direction between the initial and final wave vectors. 

For the insertion of the interactions in the u.r. case, there are different possibilities. 
The one adopted here is to insert the interaction T into the u.r. equation 

(xP+ - p + )i) » 0, (55) 

that becomes 

(p+ - xP + + P + T)tp = 0. (56) 

When this equation is written in terms of the rescaled variable £ = P + Z- we get the 
solution 

ip = exp(^-ixi + i T{b,£')d£^ (57) 
appearing in eqf52] (rotated kz — > kz ± (q/2)x). 
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Evidently, the previous u.r. and n.r. solutions are equivalent, once the interaction 
factors V and T are correctly related each other. 

The tricky point [83] is that in both cases one assumes that some observable, 
associated with the longitudinal highly relativistic motion, is constant. In the n.r. case 
it is k, in the u.r. case, it is x. "Constant" means "conserved during all the development 
of the motion". This way, the corresponding operator is substituted by its eigenvalue 
in the wave equation. 

Multiple interactions may initially convert the initial quark state with mass m into 
a more complex intermediate state with invariant mass M (e.g., a quark plus a quark- 
antiquark pair). If this state is converted again into a single quark state, it contributes to 
elastic scattering. In highly relativistic conditions we have changes of p z (in intermediate 
states) 

Ap, , M < M up (58) 

Pz 

where M up is a reasonable upper bound for the spectrum of intermediate states that 
may play some role. 

The classical Glauber regime requires M up « p z . In this case p z ~ constant in 
the interactions, the 0(p z ) term has no dynamical role, and the really relativistic part 
of the quark wavefunction disappears from the problem. 

The Gribov regime is the one where in the limit p z — > oo, M up = ap z with small 
but nonzero a. Then we have finite 0(a) changes in x during the intermediate stages 
of the process, invalidating the above assumption that x is a constant of motion. 

Since a full inclusion of inelastic channels is a difficult and unsolved problem, such 
channels are normally taken into account via two techniques: the first one is to introduce 
a discrete and limited set of intermediate non-elastic channels |8 7] . The other one is to 
rewrite the interaction terms in a non-hermitean form. 

The latter treatment is formally equivalent to a multi-channel treatment, as far as 
one is interested in the evolution of the elastic channel only. For nonrelativistic nuclear 
physics, this was demonstrated by Feshbach|82j. The scheme employed by him is quite 
abstract, and it can be generalized practically to any algebraic set of coupled equations 
describing the evolution of a state vector. So this equivalence may be considered general. 
The drawback of this generality is that it is practically impossible to use Feshbach's 
arguments to calculate the form of the effective interaction from first principles. 

So, in this work I had to guess a general form for the interaction operator T, and 
fit it on elastic proton-proton scattering according to eqEMl 

A last relevant point is that the treatment via the single channel effective potential 
does not transport the kinematics of the rescattering process out of the Glauber region 
(the price for this is that momenta are complex). Rescattering in the Glauber region 
does not create problems with factorization, at least in the Drell-Yan case for which an 
exhaustive discussion of the problem exists [72j [73] . 
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